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O ' ABSTRACT 

As increasingly precise information about the spectrum of the cosmic microwave 
background fluctuations is gathered with balloon and satellite experiments, interest has 
grown in foreground sources of opacity affecting these observations. One potentially 
important source is electron scattering produced by a post-recombination luminosity 
source, which would significantly attenuate the higher harmonics in the spectrum. If 
such an ionization source exists, then it would also heat the universe, hence increasing 
the Jeans mass and suppressing early gravitational structure formation. Here we 
consider the effects of such heating. We concentrate on one type of ionization source: 
luminosity generated by accretion onto primordial compact objects. We show that 
if such objects generate enough luminosity to affect the CMB power spectrum, then 
they would produce enough heat to prevent the formation of la collapsed objects until 
z ~ 5, significantly less than the redshift at which baryonic collapse could otherwise 
occur. Such processes would leave signatures detectable by upcoming instruments such 
as NGST, SIRTF, and SWIFT. 
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1. Introduction 

The advent of a data-rich era in cosmology is expected to produce dramatically enhanced 
understanding of the early universe, from the initial spectrum of perturbations and their processing 
before radiation decoupling, to the formation of the first galaxies and stars. The high angular 
resolution power spectra of the cosmic microwave background (CMB) gathered by current and 
future satellite, balloon, and ground-based experiments are expected to provide precise information 
about many cosmological parameters, such as the total curvature of the universe and the baryonic 
mass fraction. These power spectra are also expected to be informative about a number of 
foreground sources, such as clusters of galaxies and other large-scale structure, and about the 
epoch of ionization of the intergalactic medium (e.g., Tegmark et al. 2000). 



Early results from BOOMERanG and MAXIMA data suggested, surprisingly, that the second 
acoustic peak is significantly weaker than expected (de Bernardis et al. 2000; Hanany et al. 2000). 
More recent results from BOOMERanG (Netterfield et al. 2001), DASI (Pryke et al. 2001) and 
possibly MAXIMA (Lee et al. 2001) are consistent with the standard ACDM cosmology. However, 
the systematic uncertainties in current data are significant enough that it is still not possible to tell 
with precision whether the higher harmonics are at the expected amplitude; this will need to wait 
for instruments with absolute calibration such as MAP and Planck (e.g., Tegmark et al. 2000). 
Therefore, despite the encouraging agreement with expectations, it is worthwhile to calculate the 
impact that other effects have at both high and low redshifts. 

The CMB spectrum can strongly constrain ionization at high redshifts, but some ionizing 
processes that may have only subtle effects (but nonetheless measurable with MAP and Planck) at 
z ~ 1000 can have important consequences at lower redshifts. One such effect within the standard 
model has been suggested recently by Peebles, Seager, & Hu (2000): if there exists a radiation 
source in the universe around the time of decoupling (z ~ 1000), the resulting enhanced density 
of free electrons at high redshift would Thomson scatter the CMB and suppress the second and 
higher peaks, in addition to moving the first acoustic peak to larger angular scales. Such sources 
could have many additional cosmological effects. In particular, if scattering is posited to suppress 
the second acoustic peak (Peebles et al. 2000) then the luminosity responsible for increased 
ionization would also heat the universe. As a result, the Jeans mass at early epochs would be 
raised, which in turn could affect the redshift of collapse of the first nonlinear structures. 

One candidate for an early ionization source is accreting compact objects that exist by a 
redshift z ~ 1000, e.g., black holes formed during the quark-hadron phase transition. The first 
detailed investigation of the effects of such objects was performed by Carr (1981). He showed 
that the radiation generated by accretion of ambient gas onto primordial black holes could have a 
significant effect on the thermal and ionization history of the universe, potentially even preventing 
the universe from entering a neutral phase. Given large-amplitude baryonic fluctuations, massive 
black holes could have formed soon after recombination; Gnedin & Ostriker (1992) and Gnedin, 
Ostriker, & Rees (1995) considered the possible consequences for thermal, ionization, and 
nucleosynthesis histories of early injection of an AGN-like spectrum from such accreting sources. 

An update to Carr's analysis was performed by Miller (2000), who used the existing limits 
on electron scattering optical depth between decoupling and the current epoch to place an upper 
bound on the contribution of primordial compact objects to dark matter. This limit is based on 
the effect such scattering would have on the angular power spectrum (limits based on distortions 
of the energy spectrum are much weaker; see, e.g., Griffiths, Barbosa, & Liddle 1999). Miller 
found that the product Qco£-i{M/10 M & ) < 4 x 10~ 3 (in the notation of the current paper, 
£ < 4 x 10~ 3 ), where Qqo is the fraction of closure density in compact objects, e_i is the accretion 
efficiency divided by 10%, and M is the typical mass of the compact object. In this paper we show 
that current data imply ^<5x 10~ 4 , which corresponds to a limit on optical depth that is three 
times stronger than used by Miller (2000) based on earlier data. 



Here we explore the effect of ionizing sources on the gravitational structure formation history 
of the universe. We show that if accreting primordial compact objects generate enough ionization 
to significantly affect the second and third acoustic CMB peaks, then they raise the temperature 
of the universe to T>500 K when z ^30. As has been discussed by a number of authors (e.g., 
Ostriker & Gnedin 1996; Valageas & Silk 1999; Haiman, Abel, & Rees 2000; Gnedin 2000), an 
increase in temperature can dramatically alter the progress of structure formation. We show that 
in the present scenario, the Jeans instability threshold redshift for an n — a density fluctuation is 
reduced by ~ 5 — 7, compared to the threshold redshift when no extra heating sources are present. 
For example, 1 — a fluctuations can collapse only at z ~ 5. 

We find that the parameter £ = £lco € ~i(M / 10 Mq) governs both the distorting effect on 
the CMB power spectrum and the formation of large-scale structure. Because the two effects 
are therefore linked, we do not need to assume anything about the accretion efficiency e for the 
compact objects; instead, we may simply say that if accretion onto compact objects has a given 
effect on the CMB power spectrum, it will have a related (and much stronger) effect on the 
redshift at which nonlinear density peaks may form in baryons. If the accretion efficiency is low 
(e.g., if the flow is governed by an ADAF, CDAF, or wind-dominated solution, for which e may 
be less than 10~ 4 ; see, e.g., Narayan, Mahadevan, &; Quataert 1998; Ball, Narayan, & Quataert 
2001; Blandford & Begelman 1999), all effects are proportionally reduced, but what we emphasize 
here is the relation between possible subtle effects at high redshift and significant effects at lower 
redshift. 

In § 2 we calculate the ionization fraction produced by an accreting compact object beyond the 
immediate vicinity of the compact object ("HII region") in the surrounding ambient region which 
contains many such objects. In order to allow for nonlinear mass concentrations, we generalize the 
analysis of Miller (2000) to include accretion in a region with density differing from the average 
density of the universe. In § 3 we compute the temperature of the HII and ambient regions by 
balancing accretion heating with Compton cooling off the CMB, also assessing the importance of 
bremsstrahlung, atomic cooling, and molecular cooling. In § 4 we determine the heating implied 
if Thomson scattering is responsible for reducing the second peak to presently-observed levels. 
We then evaluate the effect this level of heating would have on the formation of self-gravitating 
structures. In §5, we discuss the effects of more general ionization and heating, and present our 
conclusions. 



2. Ionization Fraction 

In this section we focus on accreting primordial compact objects as the source of early 
ionization and heating. In § 5 we will revisit the heating effects of more general early ionization 
sources. The physical picture is of compact objects of typical mass M capturing matter via 
Bondi-Hoyle accretion from a surrounding medium. For generality we consider a medium with an 
ambient Hydrogen number density n am b and temperature T am b that may differ from the average 



density n and temperature T of the universe. This will allow us to treat accretion in underdense 
or overdense mass concentrations. The luminosity generated by accretion will produce an HII 
region immediately around each compact object, in which the ionization fraction x ~ 1. Let the 
number density and temperature inside the HII region be nnn and Thu, respectively. As discussed 
by many authors (e.g., Silk 1971; Carr 1981; Miller 2000), the hardness of the spectrum produced 
by an accreting compact object implies that its ionization effect is felt over much greater distances 
than the radius of its Stromgren sphere, and the resulting ionization fraction decreases much more 
slowly with distance than is typical around an early- type star (as a power law, not exponentially). 

Miller (2000) calculated the average ionization fraction for a differential luminosity spectrum 
typical of accreting black holes in our Galaxy, dL{E) / dE oc E^ 1 exp(— E/E max ) (the range of 
spectra observed from stellar-mass black holes all give similar answers for the ionization fraction; 
see Miller 2000 for a more detailed discussion). Typically E max ~ 100 keV. Over a region large 
compared to the radius of the Stromgren sphere, the average ionization fraction at redshifts 
z < 100 is 
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Here £"0=13.6 eV is the ground-state ionization energy of hydrogen and the number of compact 
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objects per cosmological volume is (4-7T-R 3 /3) 1 . In terms of the fraction J7 co of closure density in 



compact objects and the redshift z, 

R^ - 8 X 10- cm(l + ,)- ( J£.) "* (==*) - 1>3 itf* , (2) 

where we adopt Hq = 70km s -1 Mpc~ . The Stromgren radius is here defined as 

V 37ran amb^max; 

Here a = 1.5 x 10 _12 T3~ ' cm 3 s _1 is the recombination coefficient to the n > 2 states of hydrogen 
(Hummer 1994); recombination to the n = 1 state produces photons that are absorbed almost 
immediately, taking the standard "Case B" assumption. We scale the temperature as T = 10 3 T3 K, 
because we find typical temperatures T ~ 10 3 K at the redshifts of interest. This form for the 
recombination coefficient is accurate to better than 10% for all T > 10 2 K, and to ~ 1 — 2% for 
10 3 K< T < 10 5 K (Hummer 1994). The Bondi-Hoyle accretion rate is M = 4vrA s (GM) 2 pc7 3 , 
where for a 7 = 5/3 gas the accretion eigenvalue is X s = 0.25. This produces a luminosity per 
source L given by 
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adopting a cosmological baryon density p = 3.6 x 10~ 31 g cm _3 (l + z) 3 from Tytler et al (2000). 



The accretion efficiency is defined as L = 0.1e_iMc 2 . The Stromgren radius is then 

(5) 
Pressure balance between the ambient medium and the HII region implies nunTmi = n am bT am b. 
Collecting factors and defining £ = e^\(M /IOM & )Q CO as a parameter quantifying the generation 
of luminosity by accretion, the average ionization fraction in the ambient medium for z < 100 is 
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The factor in brackets depends only on the spectrum, and is very close to unity. We note that 
the result from equations (|l]) and (||) is a factor [8/ln(.E max /E'o)] ~ 1 times the result obtained 
by equating the total ionization at a mean rate n co L / '(3Eq) with recombinations at a mean rate 
an^ mb x 2 . Because L oc n for energy generated via cosmic accretion, the resulting mean cosmic 
ionization is density independent except indirectly through the temperature. 



3. Balance of Heating and Cooling 

The formation of baryonic structure depends in part on whether the nonlinear mass at a 
given redshift exceeds the Jeans mass at that epoch. In § 4 we will derive the nonlinear mass. 
The Jeans mass depends on the typical temperature of the regions that may undergo collapse. We 
therefore need to compute this temperature in the presence of accretion heating and ionization. 

The temperature T am b of the ambient medium is determined by the balance of heating 
by accretion with cooling by a variety of processes. Since ionization is not perfectly efficient 
(i.e., only a fraction of available energy goes into ionization), heating by accretion is inevitable. 
Cooling depends on a number of processes; candidates include Compton cooling, bremsstrahlung, 
atomic recombination, and molecular cooling. In addition to considering ionization and thermal 
balance in the ambient medium, we must also evaluate the latter in the HII region. This is 
important because, at least for heating by accretion, the luminosity per compact object depends 
on both the temperature and the density inside the HII region (the Bondi radius is smaller 
than R s for M <10 4 M Q ). We first compute the HII region and ambient temperatures Thu, 
Tumb by balancing accretion heating and Compton cooling in each region separately, assuming 
pressure balance applies at the interface. We then demonstrate self-consistency by showing that if 
ionization is dominated by radiative ionization (as opposed to collisional ionization) then cooling 
by bremsstrahlung, atomic recombination, and formation of H2 can be ignored. Only Compton 
cooling is effective. 

Balance of accretion heating and Compton cooling. - - In general, the heating rate is related 
to the ionization rate, with the heating rate per unit volume T = £tC(1 — x ) n H for <?t the typical 
excess (thermal) energy in each ionization, and C, the ionization coefficient. As discussed by Carr 



(1981), through ionization balance the heating rate may be related to the recombination rate as 
r = E^axnun e . Numerically, the heating rate at z < 100 is thus 
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where the number densities here (and everywhere else in this paper) are measured in cm~ 3 . Carr 
(1981) estimates that there is ~ 7 eV of heating per ionization within the HII region, in which 
primarily soft photons interact due to their higher absorption cross sections. In the ambient 
medium far from the source, the soft photons have already been absorbed and hence the photons 
are harder. The number of ionizations per primary of energy E is then ~ E/3Eq for typical photon 
energies ~1 keV or higher (Dalgarno, Yan, & Liu 1999), so that £ T ~ E(E/3Eq)- 1 - E = 2E . 

After decoupling but prior to structure formation the entire universe is nearly transparent 
to CMB photons. Scattering of these photons off of free electrons provides a volume cooling rate 
4a T aT£ MB k(T - TcMB)n e /{m e c), or 

Acmb = 5.6 x 10~ 33 erg cm" 3 s" 1 n e (l + z) 4 (T 3 - Tcmb, 3 ) • (8) 

At the redshifts z <C 1000 relevant to our heating calculation, the temperature in the HII region 
or ambient medium is much greater than the temperature of the CMB and hence the factor in 
parentheses is approximately T 3 . Inside the HII region, we set tih and n e to nmi, x ~ 1, and T to 
Thu; in the ambient region we set n# to n am b, n e to xn am b, and T to T am b, and evaluate x from 
equation <^). From pressure balance with the surrounding ambient medium, nHil^HJI = ^amb^amb- 

Equating the accretion heating and CMB Compton cooling rates (eqs. plS) in the HII region, 
we find 

-4/iW/li / "-amb ^ 4/11 fS^miy 11 
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Similarly, thermal balance in the ambient region yields 

T amb , 3 = 70(1 + zY^x^ (^f ) 4A (%f ) ^ • (10) 

Combining this with our expressions from equation (y) for the ionization fraction x in the ambient 
medium, and equation (0) for the temperature of the the HII region, we solve to obtain 
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The exponents are not exact, but are close approximations to the formally derived exponents. 
We neglect possible differences between £x in the HII and ambient regions. Equations (11)— (13) 
apply at redshifts where T ^> Tcmb! as z approaches decoupling, the effects of accretion become 
negligible. Figure 1 shows the high-redshift ambient ionization boost for ( = 5x 1CP 4 (which we 
find in § 4.1 gives a good fit to the second- year BOOMERanG and DASI data), compared to the 
standard model in which there is no extra ionization and hence £ = 0. 

We now consider three other cooling mechanisms: atomic recombination, bremsstrahlung, and 
molecular cooling by formation of H2 . For each of these we assume that the ionization fraction x 
is determined by the radiative ionization balance calculated earlier. In very dense or hot media, 
collisional ionization will raise the ionization fraction and therefore increase the effectiveness of 
these cooling mechanisms. 

Atomic recombination. — In solving for the ionization fraction by equating the recombination 
rate with the radiative ionization rate we have already effectively taken atomic recombination 
cooling into account. This is because the thermal energy, kT < 1 eV, is much less than the 
ionization energy of 13.6 eV. Therefore, recombination plays a minor role in cooling the medium. If 
the temperature were high enough that collisional ionization were important, then recombination 
might play a major role, but here it can be neglected. 

Bremsstrahlung. — The bremsstrahlung cooling rate is 

Abremss = 4.4 x KT^erg cm" 3 s^xum^ 2 . (14) 

The ratio of bremsstrahlung cooling to the heating rate is thus ~ 6 x 1CP 4 T 3 , which is small for 
the present situation. 

Molecular cooling. — No metals are present in the early universe, so the only ways to form 
H2 (by far the dominant molecule) are the two-stage processes (1) H + e - — > H~ + 7 followed 
by H~ + H -> H 2 + e", and (2) H + + H -> H^ + 7 followed by H% + H -> H 2 + H+. Puy et 
al. (1993) give rates for these reactions, and (1) dominates by at least two orders of magnitude 
over the temperature range of interest. Moreover, the second step of process (1) is effectively 
instantaneous compared to the first, so the reaction rate for H + e~ — ► H~ + 7 determines the 
overall rate of H2 formation. This rate is a mo i = 10 -15 s _1 cm 3 T3. Assuming as an upper bound 
that all ~10 eV released in this process escapes to infinity, the cooling rate in the (low-ionization) 
ambient medium is 

A mo i = 10 eV x a mo in e n amb . (15) 

The ratio of this molecular cooling to the mean accretion heating is 2.7 x 10~ 4 x -1 T 3 . From 
equations (|i~l"l) and (|i~3|), this ratio becomes 

A moZ /r = 0.5(1 + z)- 1 (2|*) {^f) ■ ^ 

This is generally less than unity for the redshift range of interest, for regions that are not strongly 
overdense compared to the universe as a whole. Once a region becomes significantly condensed 
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Fig. 1. — Fractional ionization a; as a function of redshift, for accretion by primordial compact 
objects (£ = 5 x 10~ 4 , the best fit to the BOOMERanG and DASI data; see § 4.1), and for no extra 
ionization source (£ = 0). Other cosmological parameters are chosen as in § 4.1. The relatively 
small difference at high redshifts implies that effects on the CMB power spectrum are small (albeit 
measurable) , but the large difference at lower redshift is one consequence of the significant heating 
effect such sources can have during the epoch of nonlinear structure formation. 



{n&mb/n S> 1), molecular cooling begin to dominate other cooling terms (Lepp & Shull 1984). For 
the purposes of assessing the onset of collapse, however, we will simply use the result expressed 
in equation (11) obtained in the approximation that accretion heating and CMB Compton cooling 



are in balance. 



4. Effects on the Formation of Structure 

With the ambient temperature T am b that we derived in the previous section, we can calculate 
the Jeans mass within a region of ambient density n am b if we know the accretion parameter 
£ = e_i(M/10MQ)Sl co . We can then compare the Jeans mass as a function of redshift to the 
nonlinear mass scales at that redshift and hence evaluate, as a function of £, whether baryonic 
condensations within nonlinear dark matter peaks can collapse. In this section we first compute 
the maximum value of £ compatible with current CMB data, then derive the nonlinear masses and 
determine the epoch of structure formation. 



4.1. Relation to the CMB power spectrum 

We use the code CMBFAST (version 4.0; Seljak & Zaldarriaga 1996 and subsequent papers) to 
generate a power spectrum and compare it to the best current signal to noise data, the second-year 
BOOMERanG (Netterfield et al. 2001) and DASI (Pryke et al. 2001) data. As an illustrative fit, 
we fixed all of the cosmological parameters except for £ and then fit for £. We assume fi& = 0.05, 
^cdm = 0.29, £1a = 0.66 (and hence assume a flat universe), Q u = 0, Hq = 66 km s _1 Mpc -1 , 
Tqmb(z = 0)=2.726 K, Yn e = 0.246, three generations of massless neutrinos and none of massive 
neutrinos. These numbers are consistent with fits of recent CMB data that include priors based 
on large scale structure, Cepheid measurements of the Hubble constant, and Type la supernovae. 
This fit is also consistent with big bang nucleosynthesis. We also assume an n = 1 initial 
perturbation spectrum and adiabatic initial conditions, and do not include tensor perturbations. 
We included the additional ionization by modifying the routine "recfast.f". At high redshift, the 
temperature of the HII region is no longer much greater than the temperature of the CMB, and 
hence for high z we solve for the ionization numerically without using this assumption. 

We plot the resulting curves for no accretion (£ = 0), £ = 5 x 10 -4 , and £ = 10~ 3 against 
the second-year BOOMERanG and DASI data in Figure 2. The vertical bars are the estimated 
statistical uncertainties in measurement. Netterfield et al. (2001) estimate that there is a 20% 
systematic uncertainty in the points, and Pryke et al. (2001) estimate that there is a 16% 
uncertainty in the power measurements, but for this figure we have plotted the raw output from 
CMBFAST against the data, including only statistical errors. Note that there are also substantial 
systematic errors contributed by beam uncertainties, particularly for £ > 400 (e.g., Netterfield 
et al. 2001), which have not been included in this figure. Clearly, the power spectrum is fairly 
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sensitive to £ in this range. To get a rough idea of the goodness of fit, we made a simplistic 
comparison of the model with the data at the center of each £ range, allowing an overall shift up 
to ±20% but not varying the data to accommodate beam uncertainties. The resulting \ 2 values 
compared against the BOOMERanG data are: 39.1 for £ = 0; 26.2 for £ = 5 x 10~ 4 ; and 65.7 for 
£ = 10~ 3 . In each case there are 17 degrees of freedom, since there are 19 data points and £ and 
the overall scale are varied. The best value of £ = 5 x 10 -4 gives a fit acceptable at the 7% level, 
not including the beam errors, which make all of the curves significantly more compatible with 
the data. We note that the best-fit value of £ = e_i$l co (M/lOM0) is compatible with a low value 
of the accretion efficiency e and/or a low value of the fraction of closure mass in compact objects 
il co , provided the compact objects are relatively massive. 



4.2. Growth of density perturbations 

When density perturbations are linear, the relative amplitude D(z) of a density perturbation 
5 = (p — p)/p evolves with the overall scale factor a = (1 + z)^ 1 of the universe according to 

m-m.-)^ £)*£*/ (it) 

(Heath 1977). Here E(z) = H(t)/H = [O m (l + zf + Q r (l + z) 4 + fi A + fi fc (l + z) 2 ] 1/2 , and 
Q m , Q r , Q\ denote the present contributions to O from matter, radiation, and the cosmological 
constant, and f2& = 1 — {Vt m + Q, r + Q\) is the overall curvature (our naming convention corresponds 
to that of, e.g., Peebles 1993). For the regime of interest, we shall neglect the £l r term. The power 
spectrum of small amplitude fluctuations at redshift z thus evolves under linear growth according 
to 

-D(zr 2 



P(k,z) 



MO) 
where k is the comoving wavenumber of a fluctuation. 



P (k) , (18) 



To obtain the power spectrum shape and normalization, we follow the development of 
Bunn & White (1997) and Hu & Sugiyama (1996). We assume that the power law slope of 
the initial perturbations is n = 1, that the universe is flat S\ = 0, and that k, the comoving 
wavenumber, is measured in Mpc -1 . Bunn & White (1997) write the present-day power spectrum 



as Ag(fc) = k 3 P {k)/27r 2 = S 2 H (jA-) T 2 (k), with 5 H ~ 2 x W~ 5 nJ- 73 . The transfer function 



,3+n 

,^o7 
is given by (Bardeen et al. 1986) 

T(q) =T(±) = ln(1 2 + 3 ' g 34g) [l + 3.89, + (16.1,) 2 + (5A6gf + (6.71g) 4 ] ^ (19) 

where hT ~ Q m h 2 . From Carroll, Press, &; Turner (1992), a good approximation to the present 
growth factor is 

m - r „ 4/7 - — , 50 : — - — , , (20) 
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Fig. 2. — Comparison of CMBFAST output, using the compact object accretion model, to the 
second-year BOOMERanG data (Netterfield et al. 2001, solid error bars) and DASI data (Pryke et 
al. 2001, dashed error bars), including statistical but not systematic errors. Here only £ is varied, 
and other cosmological parameters such as Hq and Q\ are fixed at the values indicated in the text. 
When a ±20% shift is allowed, £ = and £ = 5 x 10~ 4 fit well, but £ = 10~ 3 still gives too low 
a first peak. Note that the amplitude of the second peak is quite sensitive to £. The best value, 
£ = 5 x 10 -4 , gives an acceptable fit: for example, compared to the BOOMERanG data, x 2= 26.2 
for 17 degrees of freedom, including an overall systematic shift of 5% but not including beam errors. 
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so for n m = 0.3 and n A = 0.7, D(0) « 0.8. 

The potential for collapsed objects of a given mass to form depends first of all on whether 
perturbations of the corresponding physical scale have grown sufficiently to become nonlinear. The 
variance in the fluctuation spectrum at a given mass M (or comoving length scale R) is related to 
the power spectrum at redshift z by the following equation (e.g. Barkana & Loeb 2000): 



a 2 (M) = a 2 (R) = l o °°^k 2 P(k,z) 

we adopt a top-hat filter j\{x) = (sinx — xcosx)/x 2 . 
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A top-hat perturbation on comoving scale R collapses in an Einstein-de Sitter universe at 
the point when its overdensity as predicted by linear theory reaches a value 5 C = 1.686 (e.g. 
Peebles 1993); this critical overdensity has only a weak dependence on cosmological parameters. 
The typical redshift for a 1-er peak at mass scale M to collapse is therefore given by the implicit 
solution of a(M) = 5 C for z using equation (|2l|); the collapse redshift for a mass scale M residing 
in a 2-<T peak is given by solution of a(M) = S c /2, and so on. 

Whether the local baryons collapse into a given dark matter halo depends on the halo mass 
relative to the Jeans mass: baryons collapse with the halo only if the halo's gravitational potential 
(~ GM 2 ' 3 p i ' 3 ) exceeds the specific thermal energy of baryons (~ kT/fi). If we assume spherical 
density perturbations then the Jeans mass for a cloud of temperature T (and sound speed c s ) and 
total matter density p to t (including both baryons and cold dark matter) is 
1 / 2 \ 3 /2 

Mj = Q nptot (c£J * 7 x lo5 M ®o^/p^r 1/2 z%L*<!£ ( 22 ) 

for a cosmological abundance of helium. Here n am b is measured in cm -3 and p a mb is the ambient 
density in baryons only. If we assume that during collapse the ratio of total mass to baryonic mass 
in a mass concentration remains approximately constant at the average value (roughly 9), then 
(Ptot/Pamb) 2 ~ 1/3 and Mj k 2 x 10 5 MqTJ, 3 ra am t> • I n t ne standard model, where there is 
no extra source of heating, the minimum halo mass for collapsed baryonic objects including shell 
crossing by dark matter and other effects is (Barkana & Loeb 2000) 

M mm = 5xl0*M s — 1-^) U^ , (23) 



which is almost identical to the value M^ stand that the Jeans mass of equation ( |22|) would take 
on in the absence of cosmic heating (see eq. 41 of Barkana &; Loeb 2000). If instead heating by 
accretion is included, then from the results of the previous section we obtain 

Mj, mod « 8 x 10 10 M (l + 2)- 2 £ 7/12 (namb/n) 1/10 « 8 x 10 8 M e 7/12 [(l + ^)/10]- 2 . (24) 

Note that this mass is virtually independent of n am t,. Since £ ~ 5 x 10~ 4 fits the BOOMERanG 
data, this implies Mj )moc j w 1.0 x 10 M [(l + z)/W]~ 2 . This mass is vastly greater than it would 
be sans heating, and hence has a major effect on structure formation. 
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Figure 3 exhibits this effect. Here we show, as a function of redshift, the nonlinear mass M n \ 
for l-o", 2-<T, and 3-a peaks, along with the standard Jeans mass and the Jeans mass modified by 
accretion heating (with £ = 5 x 10 -4 ). The nominal collapse redshift for a given rarity of peak is 
significantly less than it would be without heating. 



5. Discussion and Conclusions 

We have shown that if accretion onto primordial compact objects produces enough luminosity 
to increase significantly the ionization fraction of the universe at z ~ 1000, then this luminosity 
has a dominant effect on the thermal evolution of the universe at lower redshifts. In particular, the 
Jeans mass is raised by some three orders of magnitude at z ~ 10, which would substantially delay 
baryonic structure formation. Even if the ionization produced by compact objects is insufficient 
to affect the second CMB acoustic peak strongly, the moderate dependence of the Jeans mass 
in equation (£4j) on the accretion luminosity parameter £ implies that the presence of even a 
relatively small number of (evenly-distributed) primordial compact objects could have a strong 
effect on structure formation. 

Effects on the observed cosmic background radiation other than on the power spectrum 
are very subtle. For example, the existence of extra free electrons with a temperature elevated 
above the radiation temperature will distort the photon energy spectrum. To first order this can 
be characterized by the Compton y parameter, y = J dr kAT/m e c 2 . However, even with the 
enhanced ionization and temperature, y = 2 x 10 -9 when £ = 5 x 10~ 4 , compared with the limit 
\y\ < 7 x 10~ 6 from COBE measurements (Mather et al. 1999). Thus, energy spectral distortions 
from scattering place far weaker limits on extra ionization than do limits to changes in the power 
spectrum (see also Griffiths et al. 1999). Another potential contributor to the energy spectrum is 
the luminosity from the accretion itself. However, this is also negligible. From equation @, at 
z ~ 10 3 the luminosity per source is roughly 10 36 ergs _1 e_i(M/10MQ) 2 . The total luminosity 
is L « 2 x 10 57 ergs- 1 e_i^co(M/10M o ) 2 = 2 x 10 57 erg S - 1 ^(Af/lOM ). Taking redshifts 
into account and assuming a radius of ~ 10 28 cm for the universe, the observed flux would be 
F~2x 10 -6 £ erg cm -2 s _1 . Even for £ = 1 this is a factor of ~ 10 3 less than the flux from a 
2.7 K blackbody. Taking into account that only a small fraction of the accretion radiation will be 
near the blackbody frequencies and that £ = 5 x 10~ 4 , the ratio of fluxes is ~ 10 7 . The direct 
radiation is unobservable. Thus, at high redshift the power spectrum places the strongest limits on 
accretion by primordial compact objects. However, as we have shown, at low redshifts the effects 
can be much more dramatic. 

An important difference between ionization by accreting primordial objects and ionization by, 
e.g., the first generation of stars or quasars is that the latter only exist after ordinary structure 
has already formed, whereas the former ionize and heat the universe independent of any structure 
formation. Hence, accretion onto primordial objects, or any other mechanism independent of 
structure formation at moderate redshift, can in principle delay any structure formation. Similar 
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Fig. 3. — Nonlinear mass and Jeans mass as a function of redshift, with and without accretion 
heating. The solid lines show the nonlinear mass M Q \ for la, 2a, and 3a peaks. The dotted line 
(Mj s tand) is the Jeans mass assuming no additional heating, and the long dashed line (Mj mo( j) is 
the Jeans mass modified by heating (£ = 5 x 10~ 4 ) compatible with the most recent BOOMERanG, 
MAXIMA, and DASI data. This figure shows that the effect of such heating is to delay baryonic 
structure formation drastically; typically the formation redshift of an na peak is 5-7 less than what 
it would be without the heating. 
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effects from early stars or quasars are self-limiting, and thus may not prevent - but nevertheless 
will regulate - further growth of structure (Ostriker & Gnedin 1996). 

More generally, any ionization source will also heat the universe. For example, Peebles et 
al. (2000) analyzed a simple model in which a uniform source of Ly a photons exists early in 
the universe. More hydrogen atoms are therefore in the n = 2 state, so they can be ionized more 
easily. The enhanced ionization produces more free electrons, which scatter the CMB background, 
moving the surface of last scattering to lower redshift and therefore moving the first CMB peak 
to lower £ (larger angular scales) and decreasing the amplitude of the second CMB peak. In their 
model, they assume that the rate of production of Ly a photons per unit volume is 

^ = e a n H H(t) . (25) 

Here n# is the number density of hydrogen, e a is a free parameter, and H (t) = a/a is the Hubble 
parameter at time t. Peebles et al. (2000) found that with this model, e a = 10 gave a good fit to 
the first-year BOOMERanG data. 

Let us assume that for each Ly a photon produced by the sources there is some heating 
as well, and that the energy that heats the universe is a fraction / of the energy in the Ly a 
photons themselves. At moderate and high redshifts with Q m = 0.3 and Hq = 70 km s _1 Mpc -1 , 
H m 1.2 x 10~ 18 s -1 (l + z) 3 / 2 ; the baryon density is n# = 1.5 x 10~ 7 cm~ 3 (l + z) 3 from Tytler 
et al. (2000). The energy of a Ly a photon is approximately 10 eV, so the volume heating rate 
associated with equation (^) would be 

r = fe a (WeV)n H H(t) 

= 3 x 10" 36 ergcm- 3 s- 1 /e a (l + z) 9/2 . 

From § 3, CMB cooling dominates over recombination and molecular cooling for the conditions of 
interest. If the ionization fraction is x and the temperature of the ambient medium is 10 3 T am b ) 3 K, 
the cooling rate is 

Acmb = 1 x 10- 39 erg cnrV^tl + z) 7 T ambj3 . (27) 

Equating the two, the temperature is T am \y^ = 3000/e a 2;~ 1 (l + z)~ 5 > 2 . If the hydrogen is 
completely ionized (x = 1), e a = 10 from Peebles et al (2000) implies a temperature of 1.5 x 10 4 / K 
at z = 20; reducing the ionized fraction raises the temperature. Thus, unless the generation of 
Ly a photons occurs almost free of heating (/ <C 0.01), the heating present if this source persisted 
at moderate redshift would raise the temperature dramatically above what it would be in the 
standard model, and hence increase the Jeans mass and delay structure formation. The only 
type of ionization source that could significantly affect the CMB power spectrum and yet not 
significantly affect structure formation is one with a very steep redshift dependence. If the energy 
generation rate per volume is e ex (1 + z) n with n > 7, heating effects at low redshift will be small. 
Otherwise, such heating will affect structure formation dramatically. 

The richness of data available from upcoming cosmological experiments will rapidly establish 
whether such ionizing and heating sources exist. If, for example, the first nonlinear baryonic 
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objects did not form until z ~ 15 or more recently, evidence of this should be readily apparent in 
data from many future satellites, including NGST (see, e.g., Robinson & Silk 2000) and SIRTF 
(e.g., Fazio, Eisenhardt, & Huang 1999), and in the redshifts and absorption line spectra from 
distant gamma-ray bursts observed with HETE-2 and SWIFT (Wijers et al. 1998; Lamb & 
Reichart 2000; Blain & Natarajan 2000). In addition, as more data with better calibrations are 
available from CMB experiments, the effects at higher £ will be more sensitive to the presence of 
any extra ionization (Miller 2000). In either case, the role of extra ionizing and heating sources in 
the early universe will be clarified greatly in the next five years. 

We are grateful to N. Gnedin for useful information, and to J. Ostriker, D. Lamb, and A. 
Loeb for helpful comments on the manuscript. We also thank the referee for suggestions that 
improved the clarity of the paper. This work was supported in part by NASA grant NAG 5-9756. 
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